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In view of the advent of large-scale neutrino detectors such as IceCube, the future Hyper- 
Kamiokande and the ones proposed for the Laguna project in Europe, we re-examine the deter- 
mination of the directional position of a Galactic supernova by means of its neutrinos using the 
triangulation method. We study the dependence of the pointing accuracy on the arrival time reso- 
lution of supernova neutrinos at different detectors. For a failed supernova, we expect better results 
due to the abrupt termination of the neutrino emission which allows one to measure the arrival time 
with higher precision. We found that for the time resolution of ± 2 (4) ms, the supernova can be 
located with a precision of ~ 5 (10)° on the declination and of ~ 8 (15)° on the right ascension 
angle, if we combine the observations from detectors at four different sites. 

PACS numbers: 14.60.Lm,13.15.+g,95.85.Ry 



I. INTRODUCTION 

The observation of neutrinos coming from the next 
Galactic supernova (SN) driven by gravitational core col- 
lapse (hereafter, SN implies the one caused by the gravi- 
tational collapse) is expected to provide very interesting 
information on the dynamics of the process, namely, how 
these stars explode and form black holes, see e.g., [2 [5]. 
Moreover, it may also shed light on some unknown neu- 
trino properties such as the neutrino mass ordering, see 
e.g., 014]. 

Since neutrinos can break free from the dense region of 
the star from which photons cannot escape, they will be 
the first messengers from the sky to inform us the occur- 
rence of the gravitational collapse. Indeed, it might be 
possible that the next Galactic SN cannot be located by 
optical observations due to obscuration. If so, observing 
neutrinos may be the only way to access its direction in 
the sky, apart from the possible simultaneous detection 
of gravitational waves • 

The possibility of determining the direction of a Galac- 
tic SN by merely using its neutrinos, has been investi- 
gated in the past 6-12 . Most of the authors consid- 
ered neutrino electron elastic scattering events in a water 
Chcrcnkov detector in order to determine the SN direc- 
tion H OH EH- According to Ref. [II], for a SN at 
10 kpc, the pointing accuracy is ~ 8° at 95% C.L. if the 
Super-Kamiokande detector is considered. This can be 
further improved to ~ 3° if gadolinium is added to wa- 
ter |13j . allowing to tag neutrinos from the inverse beta 
decay background. A megaton size water Cherenkov de- 
tector using this technique may be able to increase the 
pointing precision to ~ 1° [IT] . 
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On the other hand, the method of the arrival-time tri- 
angulation, previously discussed in Refs. 00, was read- 
ily dismissed due to the low precision on the arrival time 
of SN neutrinos expected mainly because the available 
detectors at that time were too small to register enough 
statistics for such a purpose. 

We are now, however, entering a new era of large- 
scale detectors with IceCube currently working in the 
South Pole [2], the proposals of Hyper-Kamiokande in 
Japan [15 and of the European detectors which will be 
built in the Pyhasalmi mine in Finland [16 . In view 
of this new trend, it is timely to revisit the usefulness 
of neutrino triangulation using big detectors in different 
continents, as suggested in Ref. [T7] . 

According to [17] , the IceCube detector can determine 
the arrival time of SN neutrinos with an uncertainty of 
± 3.5 ms at 95 % C.L.. In the case of a so-called failed 
SN, where a black hole is formed while the neutrino flux 
is still measurably high [IS], one expects the neutrino 
signal to terminate abruptly. As this sharp transition is 
expected to take place in < 0.5 ms [TJ5] the end-point of 
the neutrino spectrum can also be used for triangulation. 

Since the observation of the arrival of SN neutrinos 
by various detectors will be a valuable tool to alert as- 
tronomers about the occurrence of the star collapse [2U] 
allowing them to observe the light curve as early as possi- 
ble or the formation of a black hole (in the case of a failed 
SN), it is important to explore different approaches to re- 
construct the location of the SN as well as the failed SN 
in the sky. 



II. TRIANGULATION METHOD 

The distribution of SN in the Milky Way is expected 
to be concentrated in the Galactic disc. For the sake of 
discussion, let us consider the same SN distribution con- 
sidered in Refs. [STIES]- In Fig. [I]we show the expected 
SN distribution f(a,5), in the plane of equatorial coor- 
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FIG. 1: Expected SN probability distribution f(a,5) based on the model considered in Ref. [21J . shown in the plane of the 
equatorial coordinates a — 8 using the Hammer projection. Different contrast of the colors reflects the difference in probabilities 
as indicated in the legend. The position of the Galactic plane in the sky is indicated by the red curve. The location of the 
historical Galactic SN explosions are also shown with their type, when known, in parentheses. 



dinates a — 5 where a and 5 are, respectively, declination 
and right ascension, and f(a, 5) dadS corresponds to the 
probability to find a SN in the interval of (a, a + da) and 
(5,6 + d5) in the sky. The distribution function f(a,5) 
is normalized, as in [21] such that / da J Sf (a, 5) = 1 
with a and 5 given in radian. In this figure, we also show 
the location of the historical Galactic SN and SN1987A 
explosions. 

Let us consider two arbitrary detector sites Xi and Xj 
on the Earth and define the displacement vector as dij = 
Xi — Xj, and denote the SN direction in the sky by the 
unit vector n. Then the difference of the arrival time of 
SN neutrino signals between two detectors, Atij = ti — tj, 
is given by 



At; 



dij ■ n/c, 



(1) 



where c is the speed of light in vacuum. Here we will 
ignore the possible time delay due to the neutrino mass 
which can be estimated as 
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where D is the distance to the SN and m v is the neutrino 
mass. 

In this work, for the purpose of illustration of this 
method, we consider up to four different detector sites 
on the Earth, namely, Kamioka, South Pole, Pyhasalmi 
and ANDES (Agua Negra Deep Experiment Site) [23] 
(see also |24j). It is because four is the minimum number 
of detector positions to uniquely determine the SN loca- 
tion as we will see below. If we add more detector sites 
such as Gran Sasso and Sudbury, the results would be 
improved. We note that there is no strong dependence 
of the results as long as we select four detector locations 
which are well separated to each other. 



We note that ANDES is the first deep underground 
laboratory in the Southern Hemisphere, which could be 
constructed in the Agua Negra tunnels that will link Ar- 
gentina and Chile under the Andes, the world longest 
mountain range. In the Ref. [22] it is discussed the po- 
tential of a neutrino detector at ANDES for observations 
of SN neutrinos as well as geoneutrinos. 

In Fig. [2] we show the solution of Eq. (HI for the case 
where the SN occurs in the Galactic center, given by 
a = 17h42m27s and 5 = -28°55', for various different 
combinations of the four detector sites mentioned above. 
For definiteness, it was assumed that the SN neutrinos 
arrived at the Earth on March 20 th of 2000 at 12:00 UTC 
but it is straightforward to change this condition. 

From this plot, we can see that for a given combina- 
tion of two detector sites, the SN location can be con- 
strained, as expected, to a closed curve in the sky. It 
is also possible to see that if we have three different de- 
tector sites, we can restrict the possible SN positions to 
only two locations in the sky. For example, the curves 
for the Kamioka-South Pole and Pyhasalmi-South Pole 
combinations intersect in two locations, the true location 
of the SN as well as a fake solution. If we have detectors 
at four different sites, it is possible to eliminate the fake 
solutions and single out the true location in the sky, as 
shown in [7]. 

In practice, however, due to the finite resolution of 
the SN neutrino arrival time measurement, we can only 
establish the SN direction with limited precision. The 
accuracy of the determination of 9, the angle between the 
SN direction and the axis connecting two given detectors, 
can be roughly estimated as 



6 (cos 6) 



c 5(Atij) 



(3) 



Very roughly speaking, the accuracy of the determina- 
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FIG. 2: Possible solutions for the SN direction (a, S) consistent with a certain difference of the arrival time determined by the 
combinations of detectors located at two different sites. Here the true (input) position of the SN is assumed to be the Galactic 
center, with a — f7h42m27s and S = -28°55'. It is assumed that SN neutrinos are detected at the Earth on the vernal point 
on March 20 th , 2000 at 12:00 UTC. We consider four sites: Kamioka, Pyhasalmi, ANDES and South Pole, indicated as K, P, 
A and SP, respectively, in the legend. 



tion of the arrival time of SN neutrino signal at a given 
detector, <5£ a mvai, can be estimated as [5], 



We define our y 2 function as follows, 
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where r is the rise time of the arriving SN neutrino signal 
and N is the number of event that exists in the rising part 
of the pulse. For the case where the edge is really sharp 
or if the rising time is considered as zero, it will be 



St 



arrival 



JW 



(5) 



where r is the duration of the signal and N is the total 
number of observed events [8] . 

According to Ref. [T7] IceCube alone can reconstruct 
the signal onset of the SN neutrinos with a resolution 
of ±(1.0 — 1.7) ms at 1 u C.L., depending on the time 
interval that can be considered post bounce. For a de- 
tector like Super-Kamiokande the corresponding uncer- 
tainty could be about ±3 ms [5j, but this will depend on 
the total event rate as well as on the time structure of 
the SN neutrino signal. 

For the case of a failed SN, due to the abrupt termi- 
nation of the SN neutrino emission, expected by the for- 
mation of a black hole, this uncertainty could be further 
reduced. If the edge is really sharp (similar to the zero 
rise time case considered in [5]), this uncertainty may 
be reduced down to 0.3 ms [5], for a detector as large as 
Super-Kamiokande. 



III. COMBINED ANALYSIS 

In this section, we discuss the results of our combined 
analysis by considering observations of SN neutrinos at 
three and four different detector sites on the Earth. 
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At$> s (a ,S )-At\Y°(a,S) 
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(6) 



where At° bs (ao, 8q) is the arrival time difference of SN 
neutrinos to be observed (expected) for the input (true) 
SN location in the sky (ctQ, do) for the combination of z-th 
and j-th detector sites on the Earth whereas A£* bco (a, S) 
is the theoretically expected one for a given SN location 
(a, S). a At is the assumed time resolution. Note that by 
construction, the best fit values (a, S) obtained by our 
X 2 analysis are the solution of the Eq. (nj) for an input 
value of A^ bs . 

In Fig. [3] we show for the same input SN location at the 
Galactic center used in Fig. [2] what would be the angular 
resolution for (a, S) that would result from a combination 
of arrival time differences registered by three different de- 
tectors. For definiteness, we assume, for the combination 
of two detectors, that the arrival time difference reso- 
lution can be ± 4 ms (left panels) and ± 2 ms (right 
panels) . 

As expected, for all the combinations we considered, 
we obtained two solutions at different locations in the 
sky, the true solution and the fake one. We note the true 
and fake allowed regions are connected at 1 a C.L. for the 
cases shown in the lower four panels in Fig. [3] Though 
we have two solutions, in practice, the one that lay in 
the region of the Galactic disc has greater probability of 
being the true one. The fake solution can be eliminated 
by considering a fourth detector location as we can see 
below. 

In Fig. [4] we show the case where four different detector 
locations are considered for the same SN input location in 
the Galactic center for the time resolution of ± 4 ms (left 



Kamioka - South Pole - Pyhasalmi 

90° 

60° 




Kamioka - South Pole - Pyhasalmi 

90° 

60° 




Kamioka - South Pole - Andes 
90° 




Kamioka - South Pole - Andes 
90° 




Kamioka - Pyhasalmi - Andes 
90° 




Kamioka - Pyhasalmi - Andes 
90° 




South Pole - Pyhasalmi - Andes 
90° 




South Pole - Pyhasalmi - Andes 
90° 




FIG. 3: Cases where three detectors are considered to determine the position of the Galactic SN: Kamioka- South Pole- 
Pyhasalmi (first row), Kamioka-South Pole- ANDES (second row), Kamioka-Pyhasalmi- ANDES (third row), and South Pole- 
Pyhasalmi- ANDES (fourth row). Colors of purple, yellow and green indicate, respectively, the regions allowed at 1, 2 and 3 a 
C.L.. Here the uncertainty in the time difference measurement between two detectors is assumed to be ± 4 ms (left panels) ± 
2 ms (right panels). 



panels) and ± 2 ms (right panels) . With four detectors at 
different sites, it is possible to single out the true location 
of the SN. In Fig. [5] we show similar plots for the case 
of different SN input location, opposite to the Galactic 
center, a = 5h42m27s and S — 28°55'. From these plots 
we can conclude that the expected precision is A (a) ~ 
15(8)° and A(<5) ~ 10(5)° for the time resolution of ± 4 
ms (2 ms). 



IV. CONCLUSIONS 

The era of high-statistics neutrino detectors has 
started. IceCube is already operating in the South Pole 



and in a decade or so we expect to also have Hyper- 
Kamiokande in Japan as well as one very large neutrino 
detector in the future European underground laboratory 
in Pyhasalmi, and there is also a possibility to construct 
a new neutrino detector at ANDES. This makes the de- 
termination of the angular position of a nearby SN by 
comparing the arrival time of the first SN neutrinos at 
these different detector locations on the Earth an inter- 
esting possibility. 

The time resolution of the triangulation technique will 
be dominated by the smallest detector, since the precision 
of the reconstruction of the neutrino signal onset depends 
on the number of neutrinos registered by the detector 
(see Eqs. Q and (0)])). We have demonstrated that, in 
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FIG. 4: Allowed regions at 1, 2 and 3 a C.L. compatible with the combinations of the arrival time differences assuming four 
detector sites of Kamioka, South Pole, Pyhasalmi and ANDES for the SN at the Galactic center assuming the uncertainty in 
the time difference measurement between two detectors of ± 4 ms (left panels) and ± 2ms (right panels). 
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FIG. 5: Same as Fig. [4] but for a SN explosion that occurred at a location opposite to the Galactic center. 



general, one needs to combine the timing of four different 
detector locations in order to uniquely localize the SN 
using this method. 

Assuming a rather optimistic, but not impossible, 
the uncertainty on the arrival time difference between 
two detectors to be ~ ± (2-4) ms, we have esti- 
mated the angular resolution of the determination of 
the location of a SN that could occur in the Galac- 
tic center given by four detectors located at the South 
Pole (IceCube), Kamioka (Super-Kamiokandc or Hyper- 
Kamiokande), ANDES and Pyhasalmi (LENA, MEM- 
PHYS and GLACIER). We established that in this case 
the angular position can be known within ~ 5 (10)° in 
declination and ~ 8(15)° in right ascension for the time 
resolution of 2 (4) ms. 
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Appendix A: Description of the Equatorial 
Coordinate System 

The SN location can be given in the so-called Equa- 
torial Coordinate System by two angular coordinates, a 
known as right ascension and S known as declination. 
Right ascension measures the angular distance eastward 
along the celestial equator from the vernal equinox, its is 
analogous to terrestrial longitude. Usually right ascen- 
sion is not given in degrees but rather in sidereal hours, 
minutes and seconds. The vernal point is defined by 
where the celestial equator and the ecliptic intersect at 
00 h 00 m 00 s and longitude 0°. By definition the north ce- 
lestial pole corresponds to S = +90°, so it is analogous 
to the terrestrial latitude. 

This defines the unit vector no, which points in the 
direction of propagation of the neutrinos arriving at the 
Earth coming from the SN as 



n 



(n 0x ,noy,n 0z 
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where 



tiq x — — cos a sin 5 
n 0y = — sin a sin 5 
n 0y = — cos <5. 



(A2) 



Let us assume that a detector positioned at the i-th 
site on the Earth is localized at a certain time t by the 
following vector 



[%ii Vii Zi)j 



(A3) 



with coordinates 



Xi(t) — i?0 cos 4>i (t) sin &i 
yi(t) = R® sin (f>i(t) sinO,. 
Zi(t) = R® cos 6 h 



where R® is the radius of the Earth and 6i is the latitude 
corresponding to the position of the detector. The angle 
4>i{t) depends on time and can be given by 



fait) = <j>i(0) + cut- VLT 



where 4>i(0) is the longitude corresponding to the initial 
position of the detector, u> is the angular velocity of the 
daily rotation of the Earth and $7 is the angular velocity 
corresponding to the annual rotation of the Earth around 
the Sun. The time t refers to the moment of the day the 
SN explosion occurred (0 < t < 24 h), given in terms of 
the Coordinated Universal Time (UTC), whereas T, as- 
sumed to be common for all detectors, is the time elapsed 
after the vernal point when the detector received the SN 
neutrinos. 

So, if we have two detectors, say, one at site 1 and 
the other at site 2, we can, explicitly, write the observed 
arrival time difference as 
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- COS 2 (t)sin #2)^02: 

(A6) 



(A5) which constrains the possible values of a and 5. 
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